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Motivated by the observation of a gapless spin liquid state in ft-(BEDT-TTF) 2 Cu 2 (CN)3, we 
analyze the anisotropic triangular lattice S — 1/2 Heisenberg model with the resonating valence 
bond mean-field approximation. Paying attention to the small quasi-one-dimensional anisotropy 
of the material, we take an approach from one-dimensional (ID) chains coupled with frustrating 
zig-zag bonds. By calculating one-particle excitation spectra changing anisotropy parameter 
J' I J from the decoupled ID chains to the isotropic triangular lattice, we find almost gapless 
excitations in the wide range from the ID limit. This one-dimensionalization by frustration is 
considered to be a candidate for the mechanism of the gapless spin liquid state. 

KEYWORDS: gapless spin liquid, ft-(BEDT-TTF) 2 Cu 2 (CN)3, anisotropic triangular lattice, frustration, 
one-dimensionalization 
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Organic conductors are one of the fascinating materi- 
als which have low-dimensionality and relatively strong 
electron correlations. So far, various physical states 
have been observed and investigated intensively. 1 Among 
them, magnetism in the Mott insulating phase next to 
the unconventional superconductivity has been attract- 
ing considerable attention. This phase is observed in the 
family of ft-(BEDT-TTF) 2 X, where BEDT-TTF (ET) 
denotes bis(ethylenedithio)-tetrathiafulvalene and X rep- 
resents a monovalent anion. Similarities to that of high- 
T c cuprates are worthy of note. Another stimulating 
problem concerning magnetism is ground state proper- 
ties of geometrically frustrated spin systems such as a tri- 
angular lattice and a Kagome lattice. These two intrigu- 
ing issues meet in a material ft-(ET)2Cu2(CN)3, which 
is a Mott insulator having a nearly isotropic triangular 
lattice, and it has been in the spotlight of late. 

According to X H NMR measurements at ambient pres- 
sure, 2 ft-(ET)2Cu2(CN)3 shows no indication of long- 
range magnetic order (LRMO) down to 32 mK. This 
is 4 orders of magnitude below the exchange constant 
J ~ 250 K estimated from the temperature dependence 
of susceptibility. Recently, a similar result has been ob- 
tained by zero-field muon spin relaxation measurements, 
which have observed no LRMO down to 20mK. 3 These 
results suggest that a quantum spin liquid state is real- 
ized in the ground state. On the other hand, the static 
susceptibility remains finite down to 1.9 K, and spin- 
lattice relaxation rate 1/Ti shows power-law temperature 
dependence below 1 K. These imply that almost gapless 
spin excitation exists. This fact is a significant feature of 
the spin liquid phase observed in this material. 

Since Anderson's proposal of a resonating valence 
bond (RVB) state, 5 enormous number of studies have 
been made on the triangular lattice spin system. It is 
now a general view that the ground state of the isotropic 
triangular lattice Heisenberg model has LRMO, such as 
the 120° structure. 6 9 On the other hand, if one ne- 
glects the LRMO and assumes a disordered ground state, 
the mean-field theory of RVB state gives a spin-gap 



Table I. Anisotropy of effective transfer integrals in tt-(ET) 2 X. 13 
The definition of t and t' are not as usual (see the text). 



Anion X 


t'/t 


Cu 2 (CN) 3 


0.94 


Cu(NCS) 2 


1.19 


Cu[N(CN) 2 ]Br 


1.33 


Cu[N(CN) 2 ]Cl 


1.47 


Cu(CN)[N(CN) 2 ] 


1.47 


Ag(CN) 2 -H 2 


1.67 


Is 


1.72 



state with d cc 2_ 2/ 2+id :E2y -wave symmetry, which is called 
"d+id state". 10 12 This RVB state, describing an insu- 
lating spin system, corresponds to a projected BCS state 
at half-filling in which doubly occupied states are ex- 
cluded. Thus, the existing theories show that the ground 
state has LRMO in general, and if the magnetic order 
is destroyed in some reason, the d+id fullgap state will 
appear. If we regard the Mott insulating phase of n- 
(ET)2Cu2(CN)3 in low temperatures as an isotropic tri- 
angular lattice spin system, the results of NMR and sus- 
ceptibility measurements, which suggest neither LRMO 
nor spin gap, cannot be explained. 

In this letter, we pay attention to small anisotropy of 
k-(ET)2Cu2(CN)3 and propose a new possibility for un- 
derstanding its gapless spin liquid state. As shown in Ta- 
ble I, only ft-(ET)2Cu2(CN)3 has an opposite anisotropy 
among the family of ^-(ET^X studied in the past. Here, 
the effective transfer integrals t and t' are defined in- 
versely to the conventional way; t = corresponds to 
the square lattice, and t' = the decoupled chains. 
Therefore, ft-(ET)2Cu2(CN)3 has quasi-one-dimensional 
(Q1D) anisotropy rather than an isotropic triangular lat- 
tice. Considering that the pure ID spin system has no 
LRMO and gapless spin excitation, it is likely that this 
Q1D anisotropy is concerned with the formation of the 
gapless spin liquid state in ft-(ET)2Cu2(CN)3. 

Based on the above consideration, we study the 
Heisenberg model on an anisotropic triangular lattice, 
which is equivalent to ID chains coupled with zig-zag 
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bonds as shown in Fig. 1. The Hamiltonian is given by 



<M ; > 



<*,J> 



where <i,i'> and <i,j> represent the summation over 
intrachain and interchain nearest-neighbor pairs with an- 
tiferromagnetic coupling constant J and J', respectively 
(see Fig. 1). We investigate the anisotropy parameter 
range J' / J = 0.0-1.0, in which the model interpolates 
between the decoupled chains (J' = 0) and the isotropic 
triangular lattice (J' = J). 

In the following, we consider a projected BCS state 
defined as 



|p-BCS) =P G |BCS) 



(2) 



where Pq is the Gutzwiller projection operator which ex- 
cludes double occupancy and |BCS) is a BCS mean-field 
wave function. Since it is difficult to treat the Gutzwiller 
projection analytically, we apply an RVB mean-field ap- 
proximation to the Hamiltonian (1) and calculate the 
one-particle excitation spectra. To put it more con- 
cretely, we introduce mean fields A^- = (c^Cji) , = 

(cJt c it) = { c \i c ji) an d obtain its excitation spectrum by 
diagonalizing the mean-field Hamiltonian. This approxi- 
mation is equivalent to the "Gutzwiller approximation" 
which replaces the effect of the Gutzwiller projection op- 
erator with the statistical weight g s as 

(p-BCS^-S^p-BCS) = ^(BCS^-S^BCS). (3) 

In the simplest Gutzwiller approximation, the statisti- 
cal weight is given as g s = 4/(1 + 5) 2 where 5 is the 
density of holes, 15 and in the case of half-filling (5 = 0), 
g s = 4. Although double occupancy is no longer excluded 
from wave functions in this approximation, it is known 
in the research of high-T c superconductivity that the 
RVB mean-field (Gutzwiller) approximation gives quali- 
tatively good results. 

The spin operators Si • Sj in the Hamiltonian (1) can 
be rewritten by the fermion operators as 



Si - Si 



4i c ^ c h c ji) • ( 4 ) 
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Fig. 1. The anisotropic triangular lattice Heisenberg model with 
intrachain coupling J and interchain zig-zag coupling J' . 
ti, T2, T3 are lattice vectors. 



By introducing the mean fields, we can rewrite the 
Hamiltonian as 



Cfe (4 T c feT + c ki c k{) + (Afc4 T c f _ fei + h.c.) 



(5) 

except for constant terms. Here, and A& are given by 
€k = -3J£ Tl cos(/c • n) 



3 J' 



£ T2 cos(/c • r 2 ) + £ T3 cos(/c • r 3 ) 



(6) 



A*, = 3JA Xl cos(/c • ri) 



3 J' 



A T2 cos(/c • r 2 ) + A T3 cos(/c • r 3 ) , (7) 



where n = (1, 0), r 2 = (1/2, V5/2), r 3 = (1/2, - V5/2) 
as shown in Fig. 1, and 



(8) 



On the analogy of BCS theory, we obtain self-consistent 
equations at zero temperature 



2N ^ 

k 



Eh 



2N 4^ E k 



(9) 



with a quasiparticle excitation spectrum 



\Ai + |A fe | 



(10) 



We determine the order parameters A Ti , £ T . (i = 1,2,3) 
by solving self-consistent equations (9) numerically, and 
obtain the one-particle excitation spectrum E^. 

Firstly, we verify our method in ID limit (J / / J — 0). 
According to the exact solution, the ground state is a 
spin disordered state and the excitation spectrum is "des 
Cloizeaux-Pearson mode" with S = l. 16 In the present 
RVB mean-field theory, the one-particle excitation spec- 
trum becomes 



E k = 3J^~+\A^f\cosk x \ 



(11) 



in the ID limit. This clearly realizes gapless excitations 
at k x = ±7r/2. Note that this one-particle excitation 
describes a spin singlet breaking, i.e. 5 = 1/2 spinon 
excitation, whereas the des Cloizeaux-Pearson mode de- 
scribes S = 1 spin-wave (magnon) excitation. Thus, two- 
spinon excitations with k x = tt/2 and k x = — tt/2 form an 
S = 1 magnon with k x = 0. This means that the present 
gapless excitation spectrum obtained in the RVB mean- 
field theory is consistent with the exact des Cloizeaux- 
Pearson mode. 

Nextly, we show the results of < J' / J < 1 case, 
focusing on the following parameters 




Because of the SU(2) degeneracy at half-filling, 10 ' 15 these 
parameters are determined uniquely regardless of the de- 
generate ground states. Actually, the excitation spectrum 
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can be written as 
E 2 k = 9J 2 D 2 cos 2 k x 



+ 9 J' 2 D 2 23 



k x \^ky 



k x 



(12) 

Therefore, Di, D23 determine the dispersion relations 
along the chains (ti) and between the chains (t2,ts), 
respectively. Their J' / J dependence calculated in the 
system size L = 1200 (N = L 2 ) are plotted in Fig. 2. 
A notable feature is that D23 remains very small com- 
pared to Di, in spite of the comparatively large J' up 
to J' I J ~ 0.25. When D23 = the system is a pure ID 
chain. Indeed, when J' / J = 0, the right-hand side of the 
self-consistent equations of £ T2 , £ T3 , A X2 , A T3 become 
all equal to zero. As we show later, D23 is very small 
for J' / J < 0.25 and vanishes when J 1 / J — > 0. This in- 
dicates that there are scarcely any correlations between 
spins of different chains, and practically ID state is real- 
ized. As J' I J approaches unity, D23 gradually increases 
and becomes equal to D\. 

Finally, we show in Fig. 3 the J' / J dependence of the 
one-particle excitation spectra Ek in (12). We find that 
the structure of excitation spectra in < J' / J < 0.25 
has little difference from that of the decoupled chains 
(J 7 I J = 0.0). As a result, almost gapless excitations are 
realized in this wide parameter range. This means that 
practically ID state is realized, which is also expected 
from the behavior of D23 in Fig. 2. When J' / J exceeds 
0.25, the excitation gap gradually increases globally in 
the first Brillouin zone (1BZ). However, the shape of the 
whole spectrum is almost unchanged until the J' / J be- 
comes as large as about 0.6. Moreover, focusing on the 
lowest energy excitations (dark areas in the contour plot 
shown in Fig. 3), their locations in the 1BZ do not deviate 
from those in the ID limit (k x = ±tt/2) for J' / J < 0.8. 
Additionally, when k x = ±7r/2, the excitation spectrum 
E^ is independent of k y , i.e., Ek = 37^23 • This is be- 
cause the frustration of two interchain couplings (corre- 
sponding to the lattice vector T2 and T3) cancel the k y 
dependence. This fact is rather important, since it indi- 
cates that the excited quasiparticles along the k x = ±tt/2 
lines feel free to move along the k y direction. This is the 
same condition as in the ID limit, except for the exis- 




tence of a finite energy gap. 

Figure 4 shows the minimum gap energy in the 1BZ 
as a function of anisotropy J'/ J, changing the system 
size L. We can see the almost gapless excitations in the 
wide parameter range < J' / J < 0.25, as is already 
expected. It is quite natural that this behavior is simi- 




Fig. 2. Anisotropy dependence of D\ and D23 for L = 1200. Note 
that D23 is very small compared to Di in a wide range < 
J' I J <0.25. 



Fig. 3. Anisotropy dependence of the one-particle excitation 
spectra. Contour plots of the spectra are on the left, and sections 
along k y = line are on the right. The hexagons with broken 
lines represent 1BZ of the triangular lattice. Up to J' I J ~ 0.25, 
the spectra for each anisotropy are hardly distinguishable, and 
the one-dimensionality strongly remains for large J' / J . 
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lar to that of £>23, considering that the minimum energy 
excitations are located along k x = ±tt/2 for J' / J < 0.6. 
By plotting the same data for various system size, L, 
in a semi- log scale (Fig. 4), we can see a discontinuous 
jump for every size. We find that this critical value J'j J 
vanishes very slowly as (lnL) -1 . Thus, the discontinu- 
ity is an artifact of finite-size calculation. We also find 
that the minimum gap energy is finite when infinitesimal 
J' is introduced. Actually, we can fit the J' dependence 
as aJ' exp(— 6J/ J') 17 for J' / J < 0.6 as shown in Fig. 
4. Considering that the minimum gap energy is already 
about 3 orders of magnitude below J at J' / J ~ 0.25, it 
can be said that almost gapless excitation is realized in 
< J' / J < 0.25. This result is fairly suggestive com- 
pared with the previous series expansion 18 and linear 
spin wave 19,20 studies, all of which suggest a spin dis- 
ordered state in the parameter range J' / J < 0.25. 

From the above results, we conclude that there is a 
strong tendency to form a ID-like excitation spectrum 
for the triangular lattice spin system with anisotropy 
< J' I J < 0.6. Furthermore, even if the anisotropy 
is as large as 0.6 < J' / J < 0.8, we can still expect ID- 
like behavior for quasiparticles except for the existence of 
the excitation gap. Let us here discuss the relation to n- 
(ET)2Cu2(CN)3. The anisotropy of spin exchange inter- 
actions in this material can be estimated from J = At 2 /U 
(U being the onsite Coulomb repulsion) as J' / J ~ 0.89. 
At this anisotropy, a rather large excitation gap exists as 
shown in Fig. 4. We consider two possibilities to under- 
stand the gaplessness. One is that the small gap region in 
Fig. 4 expands to large values of J' / J by some factors not 
considered in the present model. For example, If long- 
distance exchange interactions, quantum fluctuation or 
multiple spin exchange effect 14 (higher order terms of 
the Heisenberg model) suppress not only LRMO but also 
the spin gap, we can reproduce the gapless spin liquid 
state at large J' / J . These possibilities remain as future 
problems. Another possibility is that the anisotropy J' / J 
of ft-(ET)2Cu2(CN)3 deviates from the above estimation 



10"' 



O 10' 




1 o' a ^A^.QjaA66jMiaffl^ 
0.2 



Fig. 4. (Color Online) Anisotropy dependence of the minimum 
gap energy in the 1BZ (right axis) for L=60(diamond), 120(plus), 
300 (square), 600 (cross) and 1200 (triangle). The semi-log plots of 
the same quantity are also shown (left axis). The solid line is a 
fitted exponential function aJ' exp(— bJ/J'), where a = 3.50 and 
b = 1.61. We find that the observed critical behavior is an artifact 
of finite size calculation (see the text). 



due to, for example, a finite U effect. 21 If it is in the 
range J' / J < 0.25, the excitation gap is sufficiently small 
and the susceptibility behavior (finite at 1.9 K whereas 
J ~ 250 K) can be explained. 

In summary, we analyzed an anisotropic triangular lat- 
tice Heisenberg model using RVB mean-field approxima- 
tion in order to investigate the physical origin of the gap- 
less spin liquid state observed in ft-(ET)2Cu2(CN)3. We 
payed attention to the Q1D anisotropy of this material, 
and took an approach from the ID limit. As a result of 
calculations, we found that a practically ID state with 
almost gapless excitations is realized in the wide range of 
the anisotropy parameter < J' / J < 0.25. Furthermore, 
one-dimensionality remained strongly even in J' / J > 
0.25 due to the geometrical frustration of interchain cou- 
plings. We consider this "one-dimensionalization by frus- 
tration" as a candidate for the mechanism of the gapless 
spin liquid state, although the full understanding has not 
yet been achieved. 
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